Interference effects in the Coulomb blockade regime: 
current blocking and spin preparation in symmetric nanojunctions 
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We consider nanojunctions in the single-electron tunnelling regime which, due to a high degree 
of spatial symmetry, have a degenerate many body spectrum. As a consequence, interference phe- 
nomena which cause a current blocking can occur at specific values of the bias and gate voltage. We 
present here a general formalism to give necessary and sufficient conditions for interference blockade 
also in the presence of spin polarized leads. As an example we analyze a triple quantum dot single 
electron transistor (SET). For a set-up with parallel polarized leads, we show how to selectively 
prepare the system in each of the three states of an excited spin triplet without application of any 
external magnetic field. 

PACS numbers: 73.63.Rt, 85.35.Ds, 85.35.Gv, 85.65.+h 



I. INTRODUCTION 

Single particle interference is one of the most gen- 
uine quantum mechanical effects. Since the origi- 
nal double-slit experiment^, it has been observed with 
electrons in vacuum^ and even with the more mas- 
sive C*6o molecules^. Mesoscopic rings threaded by 
a magnetic flux provided the solid-state analogous^. 
Intra-molecular interference has been recently discussed 
in molecular junctions for the case of strong^— and 
weakii~— molecule-lead coupling. What unifies these re- 
alizations of quantum interference is that the travelling 
particle has two (or more) spatially equivalent paths at 
disposal to go from one point to another of the interfer- 
ometer. 

Interference, though is hindered by decoherence. Gen- 
erally, for junctions in the strong coupling regime deco- 
herence can be neglected due to the short time of flight 
of the particle within the interferometer. In the weak 
coupling case, instead, the dwelling time is long. Usu- 
ally, the decoherence introduced by the leads dominates, 
in this regime, the picture and the dynamics essentially 
consists of sequential tunnelling events connecting the 
many-body eigenstates of the isolated system. Yet, in- 
terference is achieved whenever two energetically equiva- 
lent paths involving degenerate states contribute to the 
dynamics (see Fig. [JJ). Interference survives as far as the 
splitting between the many body levels is smaller that 
the tunnelling rate to the leads since in this limit the 
system cannot distinguish between the two paths. Thus, 
in such devices, that we called interference single electron 
transistors^ (ISET) , interference effects show up even in 
the Coulomb blockade regime. They can e.g. yield a 
selective spin blockade in an ISET coupled to ferromag- 
netic leadsi 3 -. Similar blocking effects have been found 
also in multiple quantum dot systems in do« and ao^ 
magnetic fields. 

In the present paper we develop a general theory of 
interference blockade. We give in fact an a priori algo- 
rithm for the detection of the interference blocking states 
of a generic ISET. As a concrete example, we analyze the 
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FIG. 1: Interference in a single electron transistor (SET). The 
dynamics is governed by equivalent paths in the many-body 
spectrum that involve two (or more) degenerate states. 



triple dot ISET (see Section HV)) as it represents the sim- 
plest structure exhibiting interference blockade. In par- 
ticular we concentrate on the blockade that involves an 
excited triplet state and we show how to prepare the sys- 
tem in each of the three spin states without application 
of any external magnetic field. Thus we obtain an inter- 
ference mediated control of the electron spin in quantum 
dots, a highly desirable property for spintronicsi^— and 
spin-qubit applications^—. 

The method of choice for the study of the dy- 
namics in those systems is the generalized master 
equation approach for the reduced density matrix 
(RDM), where coherences between degenerate states are 
retained^ " 13 ' 24r ~—. Such coherences give rise to preces- 
sion effects and ultimately cause interference blockade. 

The paper is organized as follows: in section [II] we in- 
troduce a generic model of ISET. In section IIIII we set 
the necessary and sufficient conditions which define the 
interference blocking states and a generic algorithm to 
detect them. In section [TV] we apply the theory to the 
triple dot molecule as archetypal example of ISET. Sec- 



tion [V] closes the paper with a summary of the results 
and conclusive remarks. 



II. GENERIC MODEL OF ISET 

Let us consider the interference single electron transis- 
tor (ISET) described by the Hamiltonian: 



states. If a blocking state is the linear combination of 
degenerate states we call it interference blocking state. 

The coupling between the system and the leads not 
only generates the tunneling dynamics described so far, 
but also contributes to an internal dynamics of the sys- 
tem that leaves unchanged its particle number. In fact 
the equation of motion for the reduced density matrix p 
of the system can be cast, to lowest non vanishing order 
in the coupling to the leads, in the for m 13 i 24 i 26 : 



H — H R 



H 



leads 



fl*, 



(l) 



1 I 

P = -T[Hsys,p\ - t[H S,P] 



£tunP- 



(3) 



where H sys represents the central system and also con- The commutator with Hsys m Eq . © represents the co 



tains the energy shift operated by a capacitivcly coupled 
gate electrode at the potential V g . The Hamiltonian _ff sys 
is invariant with respect to a set of point symmetry op- 
erations that defines the symmetry group of the device. 
This fact ensures the existence of degenerate states. In 
particular, for essentially planar structures belonging to 
the D n group, the (non-accidental) orbital degeneracy is 
at maximum twofold and can be resolved using the eigen- 
values i of the projection of the angular momentum along 
the principal axis of rotation. A generic eigenstate is then 
represented by the ket \N£aE) where N is the number 
of electrons on the system, a is the spin and E the en- 
ergy of the state. The size of the Fock space can make 
the exact diagonalization of H sys a numerical challenge 
in its own. We will not treat here this problem and con- 
centrate instead on the transport characteristics. -ffi ea ds 
describes two reservoirs of non-interacting electrons with 
a difference eVb between their electrochemical potentials. 
Finally, -f/tun accounts for the weak tunnelling coupling 
between the system and the leads, characteristic of SETs, 
and we consider the tunnelling events restricted to the 
atoms or to the dots closest to the corresponding lead: 
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where c k creates an electron with spin a and momen- 
tum k in lead a = L, R, S atJ creates an electron in the 
atom or dot closest to the lead a and t is the bare tun- 
nelling amplitude that we assume for simplicity indepen- 
dent of a, k and a. 

In the weak coupling regime the dynamics essentially 
consists of sequential tunnelling events at the source and 
drain lead that induce a flow of probability between the 
many-body eigenstates of the system. It is natural to 
define, in this picture, a blocking state as a state which 
the system can enter but from which it can not escape. 
When the system occupies a blocking state the particle 
number can not change in time and the current vanishes. 
If degenerate states participate to transport, they can 
lead to interference since, like the two arms of an elec- 
tronic interferometer, they are populated simultaneously. 
In particular, depending on the external parameters they 
can form linear superpositions which behave as blocking 



sys 

herent evolution of the system in absence of the leads. 
The operator Aun describes instead the sequential tun- 
nelling processes and is defined in terms of the transi- 
tion amplitudes between the different many-body states. 
Finally, H e g renormalizes the coherent dynamics associ- 
ated to the system Hamiltonian and is also proportional 
to the system-lead tunnelling coupling. The specific form 
of H g depends on the details of the system, yet in all 
cases it is bias and gate voltage dependent and it vanishes 
for non degenerate states. 



III. BLOCKING STATES 
A. Classification of the tunnelling processes 

For the description of the tunnelling dynamics con- 
tained in the superoperator £tun it is convenient to clas- 
sify all possible tunnelling events according to four cat- 
egories: i)Creation (Annihilation) tunnelling events that 
increase (decrease) by one the number of electrons in the 
system, ii) Source (Drain) tunnelling that involves the 
lead with the higher (lower) chemical potential, iii) | (|) 
tunnelling that involves an electron with spin up (down) 
with respect of the corresponding lead quantization axis, 
iv) Gain (Loss) tunnelling that increases (decreases) the 
energy in the system. 

Using categories i)— iii) we can efficiently organize the 
matrix elements of the system component of H tun in the 
matrices: 
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where S, D means source and drain respectively and 



t+ a = (N +l,{l',T'},E'\diJN,{£,r},E) (5) 

is a matrix in itself, defined for every creation transi- 
tion from a state with particle number N and energy E 



to one with N + 1 particles and energy E'. We indi- 
cate correspondingly in the following transitions involv- 
ing tg a and t Ecr as source-creation and drain-creation 
transitions. The compact notation {£, r} indicates all 
possible combination of the quantum numbers £ and r. 
It follows that the size of t+ a is mul(iV+l, E')xmu\(N, E) 
where the function mul(A, E) gives the degeneracy of the 
many-body energy level with N particles and energy E. 
Analogously 
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accounts for the annihilation transitions. 

The fourth category concerns energy and it is inti- 
mately related to the first and the second. Not all transi- 
tions are in fact allowed: due to the energy conservation 
and the Pauli exclusion principle holding in the fcrmionic 
leads, the energy gain (loss) of the system associated to 
a gain (loss) transition is governed by the bias voltage. 
These energy conditions are summarized in the table U 
and illustrated in Fig. [2] 



AE< 


Creation 


Annihilation 


Source 


+eV h /2 


-eVb/2 


Drain 


-eV h /2 


+eV h /2 



TABLE I: Energy conditions for tunnelling transitions be- 
tween the many-body eigenstates of the system. The quantity 
AE — Ef — Ei is the difference between the energies of the 
final and initial many-body states of the system involved in 
the transition. 
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FIG. 2: Energetically available transitions from an N particle 
level. The patterned rectangles indicate the energy range of 
energetically available source (S) and drain (D) transitions 
both to states with N + 1 and N — 1 particles. The arrows 
show examples of both allowed and forbidden transitions. 



transport and can be excluded a priori from any consid- 
eration. These are states with zero transition elements 
to all other relevant states. Within the subspace with 
iV" particles and energy E the decoupled states span the 
vector space: 



The quantity AE := Ef — Ei, is the difference between 
the energy of the final and initial state of the system 
and the approximate condition < is due to the thermal 
broadening of the Fermi distributions. For simplicity we 
set the zero of the energy at the chemical potential of 
the unbiased device and we assume an equal potential 
drop at the source and drain contact. In the table U one 
reads for example that in a source-creation tunnelling 
event the system can gain at maximum ^-^- or that in 
a source-annihilation and drain-creation transition the 
system looses at least an energy of ^rp. 

From tableUone also deduces that, from whatever ini- 
tial state, it is always possible to reach the lowest en- 
ergy state (the global minimum) via a series of energet- 
ically allowed transitions. Vice versa, not all states can 
be reached starting from the global minimum. Thus, 
the only relevant states for the transport in the station- 
ary regime are the states that can be reached from the 
global minimum via a finite number of energetically al- 
lowed transitions. 



B. Subspace of decoupled states 

In the process of detecting the blocking states we ob- 
serve first that some states do not participate to the 
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where E' is the energy of a relevant state with TV + 1 or 
N—l particles respectively. The function ker M returns 
the null space of the linear application associated to the 
matrix M. 

The decoupled space T>n,e as presented in equation 
is constructed as follows. Let us consider a generic many- 
body state \iPne) with N particles and energy E and let 
v be the vector of its components in the basis \N£tE). 
The vector T^ EE ,v has thus 4 x mul(A + 1, E') compo- 
nents and consists of all possible transition amplitudes 
from \i/jne) to all possible states with N + 1 particles 
and energy E'. Consequently hei T^ EE , contains the 
vectors v associated to states with N particles and en- 
ergy E which are decoupled from all possible states with 
N +1 particles and energy E' . Analogously holds for the 
significance of ker T~^ EE , . The intersections in ([7]) and 
the condition on E' ensure that T>n.e contains only states 
decoupled at the same time from all other states relevant 
for transport in the stationary regime. We emphasize 
that, due to the condition on the energy E' , the decou- 
pled space "Dn,e is a dynamical concept that depends 
on the applied gate and bias across the ISET. The cou- 
pled space Cn,e is the orthogonal complement of Dn,e 



in the Hilbert space with N particles and energy E. The 
blocking states belong to it. 

As a first simple application of the ideas presented so 
far, let us consider the SET at zero bias. According to 
the table U the system can only undergo loss tunnelling 
events and the global energy minimum is the only block- 
ing state, in accordance with the observation that the 
system is in equilibrium with the leads and that we mea- 
sure the energy starting from the equilibrium chemical 
potential. The potential V E of the gate electrode de- 
fines the particle number of the global minimum and, by 
sweeping V g at zero bias, one can change the number of 
electrons on the system one by one. This situation, the 
Coulomb blockade, remains unchanged until the bias is 
high enough to open a gain transition that unblocks the 
global minimum. Then, the current can flow. Depend- 
ing on the gate this first unblocking transition can be 
of the kind source-creation or drain-annihilation. Corre- 
spondingly, the current is associated to N ■<->• N + 1 or 
N -H> N — 1 oscillations, where N is the particle number 
of the global minimum. 



C. Blocking conditions 

At finite bias the condition which defines a blocking 
state becomes more elaborate: 

1. The blocking state must be achievable from the 
global minimum with a finite number of allowed 
transitions. 

2. All matrix elements corresponding to energeti- 
cally allowed transitions outgoing from the block- 
ing state should vanish: in particular all matrix 



elements corresponding to Ej — E\ 
for \E t -E h 



^block 
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^k only the ones corresponding 



'Jf — ^blocki -s 2 

to the drain-annihilation and source-creation tran- 
sitions. 

The first condition ensures the blocking state to be 
populated in the stationary regime. The second is a 
modification of the generic definition of blocking state 
restricted to energetically allowed transitions and it can 
be written in terms of the tunnelling matrices T^ EE , 

and Tft EE i- For each many-body energy level \NE), the 
space spanned by the blocking states reads then: 



with 
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In Eq. |H]we introduced the matrices Tp = (0, 1) T and 
Ts = (1,0) T with 1 being the identity matrix and 
the zero matrix, both of dimension 2 x mul(A + 1,25') 
for T D and 2 x mu\(N - 1,25") for T s . The energies 
E' and E" satisfy the inequalities \E' — E\ < ^p and 
E" — E < — ^Tp 1 , respectively, and Vne is the projection 
on the N particle space with energy E. 

The first kernel in B N E together with the projector 
Vne gives all linear combinations of N particle degener- 
ate states which have a finite creation transition involving 
the drain but not the source lead. This condition can in 
fact be expressed as a non-homogeneous linear equation 
for the vector v of the components in the many body 
basis of the generic TV particle state with energy E: 



1 N,EE' V ~ "' 



(10) 



where b is a generic vector of length 4 x mul(iV + 1, E') 
whose first 2 x mul(./V + 1,25') components (the source 
transition amplitudes) are set to zero. Due to the form 
of b, it is convenient to transform Eq. (|10[) into an 
homogeneous equation for a larger space of dimension 
mul(iV, E) + 2 x mul(7V + 1, E') which also contains the 
non-zero elements of b and finally project the solutions of 
this equation on the original space. With this procedure 
we can identify the space of the solutions of (jTDJ) with: 



V = V NE ker (T+ EE ,,T D ) 



(11) 



,(i) 



The second kernel in B N E takes care of the annihi- 
lation transitions in a similar way. Notice that V also 
contains vectors that are decoupled at both leads. This 
redundance is cured in @ by the intersection with the 
coupled space Cne- 

The conditions ^ are the generalization of the condi- 
tions over the tunnelling amplitudes that we gave in [l2J . 
That very simple condition captures the essence of the 
effect, but it is only valid under certain conditions: the 
spin channels should be independent, the relevant energy 
levels only two and the transition has to be between a non 
degenerate and a doubly degenerate level. Equation (0), 
on the contrary, is completely general. In appendix IA1 we 
give an explicit derivation of the equivalence of the two 
approaches in the simple case. 

For most particle numbers N and energies 25, and suf- 
ficiently high bias, Bn,e is empty. Yet, blocking states 
exist and the dimension of Bn,e can even be larger than 
one as we have already proven for the benzene and the 
triple dot ISETs^. Moreover, it is most probable to find 
interference blocking states among ground states due to 
the small number of intersections appearing in (j9j) in this 
situation. Nevertheless also excited states can block the 
current as we will show in the next section. 

The case of spin polarized leads is already included 
in the formalism both in the parallel and non parallel 
configuration. In the parallel case one quantization axis is 



naturally defined on the all structure and a in equations 
([5]) and © is defined along this axis. In the case of non 
parallel polarized leads instead it is enough to consider 
S a(T and d aa in equations (|5|) and ([5]) , respectively, with a 
along the quantization axis of the lead a. It is interesting 
to note that in that case, no blocking states can be found 
unless the polarization of one of the leads is P = 1 . The 
spin channel can in fact be closed only one at the time 
via linear combination of different spin states. 

A last comment on the definition of the blocking con- 
ditions is necessary. A blocking state is a stationary so- 
lution of the equation © since by definition it does not 
evolve in time. The density matrix associated to one of 
the blocking states discussed so far i) commutes with the 
system Hamiltonian since it is a state with given particle 
number and energy; ii) it is the solution of the equation 
AunP = since the probability of tunnelling out from 
a blocking state vanishes, independent of the final state. 
Nevertheless, a third condition is needed to satisfy the 
condition of stationarity: 

3. The density matrix pbiock associated to the blocking 
state should commute with the effective Hamilto- 
nian H c f{ which renormalizes the coherent dynam- 
ics of the system to the lowest non vanishing order 
in the coupling to the leads: 



[Pblock, H c ff] — 0. 



(12) 



The specific form of H c s varies with the details of the 
system. Yet its generic bias and gate voltage dependence 
implies that, if present, the current blocking occurs only 
at specific values of the bias for each gate voltage. Fur- 
ther, if an energy level has multiple blocking states and 
the effective Hamiltonian distinguishes between them, se- 
lective current blocking, and correspondingly all electri- 
cal preparation of the system in one specific degenerate 
state, can be achieved. In particular, for spin polarized 
leads, the system can be prepared in a particular spin 
state without the application of any external magnetic 
field as we will show explicitly in section HV CI 



IV. THE TRIPLE DOT ISET 

The triple dot SET has been recently in the focus 
of intense theoretica l 14 i 15 ' 37 ~ — and experimental^ - — in- 
vestigation due to its capability of combining incoher- 
ent transport characteristics and signatures of molecular 
coherence. The triple dot ISET that we consider here 
(Fig. [3]) is the simplest structure with symmetry pro- 
tected orbital degeneracy exhibiting interference block- 
ade. Despite its relative simplicity this system displays 
different kinds of current blocking and it represents for 
this reason a suitable playground for the ideas presented 
so far. 



A. The model 

We describe the system with an Hamiltonian in the 
extended Hubbard form: 



+[/^(n it -i)(n^-i) 

i 

+v ^ ( n n + n u _ i) ("i+iT + n i+n - 1) i 



(13) 



ft 



where d ia creates an electron of spin a in the ground 
state of the quantum dot i. Here i = 1, . . . , 3 runs over 
the three quantum dots of the system and we impose the 
periodic condition d^ a = d\ a . Moreover n^ = d ia d i(J . 
The effect of the gate is included as a renormalization 
of the on-site energy £ = £o — eV g where V g is the gate 
voltage. We measure the energies in units of the modulus 
of the (negative) hopping integral b. The parameters that 
we use are £ = 0, U = 5 |6|, V = 2\b\. 




FIG. 3: Schematic representation of a triple dot interference 
single electron transistor (ISET). 

The number of electrons considered for the triple dot 
structure goes from to 6. Thus the entire Fock space 
of the system contains 4 3 = 64 states. By exact diago- 
nalization we obtain the many body-cigenstates and the 
corresponding eigenvalues that we present in Fig. 2] for a 
gate voltage of V g = 4.8 6/ e. In the table |H] we also give 
the degeneracies of all levels relevant for the blocking 
states analysis which will follow. We distinguish between 
spin and orbital degeneracy since the latter is the most 



important for the identification of the blocking states. 
The total degeneracy of a level is simply the product of 
the two. 



Many-body 


Orbital 


Spin 


energy level 


degeneracy 


degeneracy 





1 


1 


lo 


1 


2 


2o 


1 


1 


2i 


2 


3 


3o 


2 


2 


4o 


1 


3 


5o 


2 


2 


6 


1 


1 



TABLE II: Degeneracy of the triple dot system energy levels 
as it follows from the underlying D3 symmetry. A level Ni is 
the ith excited level with N particles. The total degeneracy 
of the level is the product of its orbital and spin degeneracies. 
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FIG. 4: Spectrum of the triple dot system for the specific gate 
voltage eV g = 4.86 chosen to favor a configuration with two 
electrons. The other parameters in the system are U = 5|6| 
and V = 2\b\, where 6 is the hopping integral between the 
different dots. 



B. Excited state blocking 

In Fig. [5] we show the stationary current through the 
triple dot ISET as a function of bias and gate volt- 
age. At low bias the current vanishes almost everywhere 
due to Coulomb blockade. The particle number is fixed 
within each Coulomb diamond by the gate voltage and 
the zero particle diamond is the first to the right. The 
zero current lines running parallel to the borders of the 
6, 4 and 2 particle diamonds are instead signatures of 



ground state interference that involves an orbitally non- 
degenerate ground state (with 2, 4, and 6 particle) and an 
orbitally double-degenerate one (with 3 and 5 particles). 
In appendix [X] we illustrate how to obtain an expression 
for the blocking states in these cases. 

The striking feature in Fig. [5] is the black area of cur- 
rent blocking sticking out of the right side of the two par- 
ticles Coulomb diamond. It is the fingerprint of the oc- 
cupation of an excited interference blocking state. Fig. [5] 
is a zoom of the current plot in the vicinity of this ex- 
cited state blocking. The dashed lines indicate at which 
bias and gate voltage a specific transition is energetically 
allowed, with the notation Ni labelling the ith excited 
many-body level with N particles. These lines are phys- 
ically recognizable as abrupt changes in the current and 
run all parallel to two fundamental directions determined 
by the ground state transitions. For positive bias, posi- 
tive (negative) slope lines indicates the bias threshold for 
the opening of source-creation (drain-annihilation) tran- 
sitions. The higher the bias the more transitions are 
open, the higher, in general, the current. 
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FIG. 5: Stationary current for the triple dot ISET. Coulomb 
blockade diamonds are visible at low biases. Ground state 
and excited state interference blockades are also highlighted. 
The temperature is kbT — 0.002|b|. The other parameters are 
the same as the ones in Fig. [3] 

The anomalous blockade region is delimited on three 
sides by transitions lines associated to the first excited 
two particle level 2±. Our group theoretical analysis 
shows that the two particle first excited state is a twofold 
orbitally degenerate spin triplet (see Table Ql]) . In other 
terms we can classify its six states with the notation 
\2i,£, S) with I = ±h being the projection of the an- 
gular momentum along the main rotation axis, perpen- 
dicular to the plane of the triple dot, and S z = —h,0,h 
the component of the spin along a generic quantization 
axis. The lo energy level is instead twice spin degener- 
ate and invariant under the symmetry operations of the 
point group D3. 

In order to identify the 2 particle blocking states we 
perform the analysis presented in the previous section 
for the 2i energy level with the gate and bias in the 
blocking region. Firstly, we find that the 2i energy level 
can be reached from 2o via the drain-annihilation tran- 
sition 2o — > lo followed by the source-creation transition 
lo — > 2i. Secondly, the space of the decoupled states 
T>2 1 is empty and the only energetically allowed outgoing 




0.4 



0.3 



0.2 



0.11 



fe 



4 4.5 5 5.5 

U|b|] 

FIG. 6: Blow up of the stationary current through the triple 
dot ISET around the 2 to 1 particle degeneracy point. The 
black area sticking out of the 2 particles Coulomb diamond 
denotes the excited states blocking. 



transition is the drain-annihilation 2i — > 1q transition. 
Thus the blocking space is given by the expression: 
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(14) 



and has dimension three. For clearness we give in the 
appendix [B] the explicit expression of T 2 _ 2 1 and the cor- 
responding vectors that span Bi x ■ Essentially, there is a 
blocking state for each of the three projection of the spin 
$2 . This result is natural since, for unpolarized or paral- 
lel polarized leads, coherences between states of different 
spin projection along the common lead quantization axis 
do not survive in the stationary limit. 

Outside the blocking region either the first or the sec- 
ond blocking state conditions arc violated. In particu- 
lar, below the lower right border the state 2i can not 
be reached from the global minimum since the lo — > 2i 
source-creation transition is forbidden while above the 
upper left (right) borders the state 2i can be depop- 
ulated towards the 3o (lo) states via a source-creation 
(drain-annihilation) transition. 



C. Spin polarized transport 

The orbital interference blocking presented in the pre- 
vious section acquires a spin dependence in presence of 
polarized leads. The lead polarization P a with a = L,R 
is defined by means of the density of states D aa at the 
Fermi energy for the different spin states: 



the system via the spin dependent bare tunnelling rates 
Tqo- = WtfD aG that enter the definition of the tun- 
nelling component of the Liouvillian £tun and the renor- 
malization frequencies u} a s z ■ We assume the leads to be 
parallel polarized so that no spin torque is active in the 
device and we can exclude the spin accumulation associ- 
ated to thati^&. 

In Fig. [7]we show the current in the excited state block- 
ing region as a function of the bias and of the (parallel) 
lead polarization P. For non-polarized leads the current 
is blocked at a single bias, while for finite values of P the 
blocking is threefold. For the same bias and polarization 
ranges we present in Fig. [5] the z component of the spin 
for the triple dot. The spin projection S z assumes, ex- 
actly in correspondence of the current blocking, the val- 
ues S z = —h,0,H, respectively, as the bias is increased. 
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FIG. 7: Current as a function of the polarization and of 
the bias voltage in proximity of the excited state interference 
blocking. The white dashed lines are defined by the conditions 
u)ss z = for S z = —h, 0, h from left to right, respectively. 




P a = 



Pat ~ Daj 



(15) 



and is taken equal for the two leads P = Pl = Pr- 
Finally, the spin polarization influences the dynamics of 



FIG. 8: Spin projection S z as a function of the polarization 
and of the bias in proximity of the excited state interference 
blocking. Notice by comparison with Fig. [7| the correspon- 
dence between current blocking and spin preparation. 

The explanation of this effect relies on the third block- 



ing condition, Eq. (JT2J , and concerns the form of the 
effective Hamiltonian introduced in Eq. ([3]). The latter 
can be written, due to the rotational symmetry of the 
system, in the form: 



H e ff = 2_^ UJ aS z L a , 



(16) 



aS z 



where L a is the projection of the angular momentum 
in the direction of the lead a and it does not depend 
on the spin degree of freedom S z . Moreover, w a s. is 
the frequency renormalization given to the states of spin 
projection S z by their coupling to the a lead. In the 
appendix [C] we give an explicit expression for ui a s, and 
L a . In Fig. |H]we plot instead ujls^ as a function of the 
bias for different polarizations. The gate is fixed at V g — 
4.8b/ e. 
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FIG. 9: Renormalization frequencies as a function of the bias 
for different polarizations P = 0,0.25,0.5,0.75 in the leads. 
The gate is V g — 4.8 b/e. At P = all the renormalization 
frequencies coincide (full line). The S z — h( — h) frequency in- 
creases (decreases) monotonously with the polarization. The 
S z — frequency is instead independent of it. 

Since the two particle ground state is totally symmet- 
ric (Ai symmetry), a three particle blocking state must 
be antisymmetric with respect to the vertical plane that 
intersects the center of the system and the drain dot. 
For this reason a blocking state is also an eigenstate of 
the projection Ld of the angular momentum in the di- 
rection of the drain lead. Consequently, the last blocking 
condition is satisfied only if: 



uss z = 



(17) 



and the effective Hamiltonian is proportional to Ld- 

For zero polarization in the leads the condition (fTT)) 
holds at the same bias for the three spin projections S z 
and the blocking state is a statistical mixture of the three 
spin projections. For polarized leads, instead, each spin 
projection is blocked at a specific bias and the spin on the 



system is controlled simply by changing the bias across 
the device. The dashed lines in Figs. [7] and \E\ represent 
the solutions of equations (fl7|) for S z = —ft, 0, ft from left 
to right, respectively. Clearly they also indicate in Fig. [7] 
the zeros of the current and in Fig. [5] the fully populated 
spin states. 



V. CONCLUSIONS 

In this paper we addressed the interference effects that 
characterize the transport through a symmetric single 
electron transistor. In particular we gave the generic con- 
ditions for interference blockade and an algorithm for the 
identification of the interference blocking states as linear 
combination of degenerate many-body eigenstates of the 
system. 

As an application of the theory we studied the triple 
dot ISET. Despite its relative simplicity, this system ex- 
hibits different types of interference blocking and it rep- 
resents an interesting playground of the general theory. 
Specifically, we concentrated on the interference block- 
ade that involves an excited triplet state. In presence of 
polarized leads we exploited the interference blocking in 
order to access each of the triplet states by all electrical 
means. 

The theory is sufficiently general to be applied to any 
device consisting of a system with degenerate many-body 
spectrum weakly coupled to metallic leads e.g. molecular 
junctions, graphene or carbon nanotube quantum dots, 
artificial molecules. In particular, the algebraic formu- 
lation of the blocking condition in terms of kernels of 
the tunnelling matrices T^, Eq. ([§]), allows a straightfor- 
ward numerical implementation and makes the algorithm 
directly applicable to complex junctions with highly de- 
generate spectrum. 
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Appendix A 

We derive here the equation (1) in [12| as a specific 
example of the general theory presented in the paper. 
That equation represents the interference blocking con- 
dition for the simplest possible configuration involving 
only a non degenerate and a doubly degenerate state. 

Let us consider for simplicity a spinless^ system and a 
gate and bias condition that restricts the set of relevant 
many-body states to three: one with A^ particles and 
two (degenerate) with N + 1 or N — 1 particles. The 
interference blocking state, if it exists, belongs to the 
N ± 1 level. There is only one interesting tunnelling 



matrix to be analyzed, namely T^ ±1 . Let us take for 
it the generic form: 



A N±1 



75i 7S2 

7-Dl 7C2 



(Al) 



where S and D indicate source and drain respectively 
and 1 and 2 label the two degenerate states with N ± 1 
particles. Js(D)i are the elements of the tg, D ^ matrices 
introduced in Eqs. §5§ and ©• 
The decoupled space reads: 



V 



N±l 



kcrT N±1 . 



(A2) 



Since the N ± 1 particles relevant Hilbert space has 
dimension 2 the only possibility to find a blocking state 
is that T>n±i = 0- In other terms: 



detT* ±1 



7S17-D2 — 7-D17S2 7^ 



(A3) 



This condition is identical to Eq. (1) ir>i2. The blocking 
state can finally be calculated as: 



Bn+i = Pjv+ikcr 



or 






JV-l 



pAT-iker 



751 7S2 1 

7-Dl 7-D2 

751 7S2 

7-Dl 1D2 1 



nc 



N+l 



(A4) 



nc 



JV-l) 



where the C/v±i is, in the relevant case, the entire space 
and the projector Vn±i simply removes the last com- 
ponent of the vector that defines the one dimensional 

kernel. 



J 



Appendix B 

We give here explicitly the T 2 _ 2i , matrix necessary 
for the calculation of the triplet blocking states and the 
associated blocking states. The states in the lo doublet 
and in the two times orbitally degenerate triplet 2i are 
labelled and ordered as follows: 





' \2 1 ,t = +h,s z = +h) 




\2 1 ,£ = +h,S z = 0) 


o,/ = 0,t) 2 


\2 1 ,i = +n,s z = -h) 


0,* = <U) ' 


\2 1 ,£ = -h,S z = +H) ' 




\2 1 ,£ = -h,S z =Q) 




, \2 1 ,£ = -h,S z = -h) 


(Bl) 


ents of the t aa matrices that compose T 2 2 i 



have thus the general form: 

t- a (s z , s' z ,e) = (i Q ,£' = o, s' z \d aa \2 } , e,s z )._ 

By orbital and spin symmetry arguments it is possible to 
show that 

t- a {S z ,S' z ,l) =tei e ^ 6 Si , S ;-* (VaJsi,f. + $su) 



where 



t = (l ,£' = 0,i\d Mf \2 1 ,£ = l,S z = 0) 



The subscript M labels a reference dot and 4> a is the 
angle of the rotation that brings the dot a on the dot M. 
The explicit form of T 2 ~ 2 1 reads: 



T, 



2,2, In 



v / 2e +l27r / 3 



/ V2e~ i27r / 3 

e -^/3 

e -^/3 

VZe- 12 "/ 3 V2e +l2 "/ 3 

V2e+ l27r / 3 V2e~ l27r / 3 





+ i2TT/3 

3 +»27r/3 



e +i27r / 3 

e +i^/3 o 

V2e +l27r / 3 



r 



e -i2ir/3 o 

e -i2ir/3 o 

V2e-^/ 3 J 



(B2) 



The rank of this matrix is 6 since all columns are inde- 
pendent. Thus C2,2i coincides with the full Hilbert space 
of the first excited 2 electron energy level. The blocking 
space jB2,2i,i reads: 



#2,2!, in = ^2 1 ker(T 22i lo ,T s ) 
where T$ reads 



(B3) 



To 



/iooooooo\ 

10 
10 

\0 1 0/ 



(B4) 
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in accordance to its general definition given in Eq. ©, 
and the projector V 2l removes the last four components 
from the vectors that span ker(T 2 ~ 2 x ,Ts). It is then 
straightforward to calculate the vectors that span the 
blocking space B2,2 1 ,i '- 



f*=2-\ 



Vl 



s/2 










v 2 



e ! e 






V o J 



V3 



V2 






V o / V o / \ e -^tJ 

(B5) 
The vectors vi , v 2 and V3 are the components of the 
blocking states written in the 2i basis set presented in 
(|B1[) . Thus, the three blocking states correspond to 
the three different projectors of the total spin S z = 
H, 0, and — h, respectively. 



where Vne = Tl mT \NmTE)(NmTE\ is the projector on 
the TV-particle level with energy E and d M(J destroys an 
electron of spin a in the middle dot M. We defined 



the function p a (x) 



-Reip 



&-(x- 



Hot, 



where 



/3 = l/k^T, T is the temperature and tp is the digamma 
function. Moreover T Q a(J , = ^-\t\ 2 D arT > is the bare tun- 
nelling rate to the lead a of an electron of spin a', where 
t is the tunnelling amplitude and D aa > is the density of 
states for electrons of spin a' in the lead a at the cor- 
responding chemical potential fi a . Due to the particular 
choice of the arbitrary phase of the 2 particle excited 
states, Was, does not depend on the orbital quantum 
number I. It depends instead on the bias and gate volt- 
age through the energy of the 1, 2i and 3 particle states. 

In the Hilbert space generated by the two-fold orbitally 
degenerate \2i£S z ) the operator L a reads: 



Appendix C 

We present here explicitly the rcnormalization fre- 
quency uJaS- and the projection of the angular momen- 
tum L a which appear in the expression of the effective 
Hamiltonian (|16|) . The frequency uj a s z is defined for the 
degenerate two particle excited level 2i in terms of tran- 
sition amplitudes to all the states of neighbor particle 
numbers: 



w aS s = - 2, 1"° a' 



a'E 



{2 1 tS z \d M(J ,V 3E d\ Ia ,\2 1 ~iS z ) Pa (E-E 2l 
(2 1 £S z \dl a ,V lE d Ma ,\2 1 -£S z )p a (E 2l -E) 



L a — 



J2<p a 



-i2<j> a 



(C2) 



where (j) a = ±^ is the angle of which we have to rotate 
the triple dot system to bring the middle dot M into the 
position of the contact dot a. For a derivation of (|C2p 
see the supplementary material of [13[ . For all degenerate 
subspaces, if no accidental degeneracy is present (like for 
our parameter choice), the effective Hamiltonian has the 
form given in ((16)) . (|C1|) . (|C2I) . with the renormalization 
frequencies calculated using the appropriate energies and 
matrix elements. 
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